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Abstract 

We study heat conduction in a harmonic crystal whose bulk dynamics is supplemented by ran- 
dom reversals (flips) of the velocity of each particle at a rate A. The system is maintained in a 
nonequilibrium stationary state(NESS) by contacts with white noise Langevin reservoirs at differ- 
ent temperatures. We show that the one-body and pair correlations in this system are the same 
(after an appropriate mapping of parameters) as those obtained for a model with self-consistent 
reservoirs. This is true both for the case of equal and random(quenched) masses. While the heat 
conductivity in the NESS of the ordered system is known explicitly, much less is known about the 
random mass case. Here we investigate the random system, with velocity flips. We improve the 
bounds on the Green-Kubo conductivity obtained by Bernardin [5] . The conductivity of the ID 
system is then studied both numerically and analytically. This sheds some light on the effect of 
noise on the transport properties of systems with localized states caused by quenched disorder. 
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I. INTRODUCTION 



We consider heat transport in mass disordered harmonic lattices with stochastic bulk 
dynamics. For the ID disordered harmonic lattice without stochasticity the effect of local- 
ization due to disorder leads, in the presence of pinning, to an exponential decay of the heat 
current as a function of the length [U [2] . In the absence of pinning the conductivity depends 
on the boundary conditions either growing as y/N or decaying as 1/ \/N j3j H] . The situation 
is very different when one adds stochasticity to the dynamics. Bernardin |5] obtained finite 
positive upper and lower bounds on the Green-Kubo conductivity of a harmonic lattice with 
periodic boundary conditions subjected to stochastic dynamics which conserves energy but 
not momentum. 

In this paper we study the heat flux of such a disordered harmonic system, both pinned 
and unpinned, connected to white noise Langevin type heat reservoirs with different tem- 
peratures at the two ends. The stochastic part of the bulk dynamics consists of random 
reversals of each particle's velocity at a rate A. The analytical as well as accurate numerical 
tractability of this model makes it a good test system to address the problem of the effect of 
noisy dynamics on transport properties of disordered systems. To the extent that stochastic 
dynamics affect phonon-phonon interactions in some rough sense similar to anharmonicity, 
this may also teach us something about the effects of the latter. The effect of interactions 
on localization and transport in disordered systems has generated much interest and has 
been studied both for the case of electron and phonon transport [21 15HT3] . A study by two 
of us [2] of the heat conductivity of the disordered ID pinned lattice with anharmonicity 
found that a small amount of anharmonicity was sufficient to cause a transition to a diffusive 
regime with a finite value of k. The question of whether the transition from an insulator 
to conductor occurs at zero or some finite small value of anharmonicity remains an open 
problem. For the noisy dynamics on the other hand [3] shows that an arbitrarily small noise 
leads to normal transport. The work of [13] suggests a transition at zero nonlinearity for 
classical systems (the quantum situation may be different). 

There have been some recent studies [HI [15] on the NESS of ordered harmonic chains 
with noisy dynamics conserving both energy and momentum, and with white noise Langevin 
baths at different temperatures at the two ends. It was shown in these studies that the time- 
evolution of the pair-correlations formed closed sets of equations. This closure is true also for 
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the energy conserving model that we study. In addition we show a mapping of the steady 
state equations for one-body and pair-correlation functions to that of the self-consistent 
reservoirs model first introduced by Bolsterli etal [TBI E] arid recently solved exactly for the 
ordered case by Bonetto etal [18] . 

The plan of the paper is as follows. In sec [IT] we define the model precisely and show 



the mapping to the model with self-consistent reservoirs. In sec III we use the method 
of Bernardin [5J to obtain improved lower and upper bounds for the Green-Kubo thermal 



conductivity of the random mass system. In sec IV we present results from numerical 
calculations as well as nonequilibrium simulations for the dependence of the heat flux in the 
random mass case on system size and on noise strength. Both the pinned and unpinned 
system are studied. For the pinned case conductivity decreases with the strength of the 
binding and increases with strength of the noise. For the unpinned case we study the effect 
of boundary conditions (BCs) [23]. Finally we conclude with a discussion of the nature of the 
NESS for this model. While it is easily checked that the NESS is not strictly Gaussian, we 
find that the one particle momentum and position distributions are very close to Gaussian 
distribution. Furthermore their values at different sites or of momentum and position at the 
same site are essentially uncorrelated. 



II. NOISY DYNAMICS AND SELF-CONSISTENT RESERVOIRS 



The results in this section apply in any dimension. For simplicity of notation we consider 
here explicitly the one dimensional harmonic chain with the Hamiltonian: 

"= E 

l=l,N 

= - jpM _1 p + q&q , 

where {qi,pi} denote the position and momenta of the particles. We have used the notation 
P = (P11P2, ■■-,Pn), 1 — (qi,Q2, •••■>Qn) an d M and $ are N x N matrices corresponding to 
masses and forces respectively. When k a > we have the pinned case and set k' = k. In the 
unpinned case, k = 0, we consider fixed, k' > 0, and free, k' = 0, boundary conditions. 

The system's evolution has a deterministic part described by the Hamiltonian above and 
a stochastic part consisting of two different processes: (i) every particle is subjected to a 
noise which flips its momentum, i.e. for the I th particle the transition pi — > —pi occurs with 
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a rate A, (ii) the particles at the boundaries I = 1 and / = iV are attached to heat baths 
with Langevin dynamics at temperatures Tl and Tr respectively. Thus the end particles 
have additional terms in their equation of motion of the form —^p a /m a + (2 , yT a ) 1 ^ 2 r] a (t), for 
a = 1,N, with {rj a (t)r] a '(t')) = 5 ata > S(t—t'), 7 is the friction constant and T\ = Tl, Tn = Tr 
are the bath temperatures. 

The master equation describing the time evolution of the full phase space probability 
density is therefore given by: 

dP(x) sr^ ~ d „ di m d 2 P , v^r„, \ ^, 

where x = (gi, q 2 , ■ ■ ■ , Qn,Pi,P2, ■ ■ ■ ,Pn) = (^1, x 2 , • • • , %2n) and d and d are 2N x 2iV ma- 
trices given by: 

' — M -1 



$ Af-T 



lr-i 
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Here T and f are diagonal matrices with diagonal elements given by Tu = Tl5i,i + TrS^n 
and fa = 7(5;, 1 + S^n) respectively. Similar to the case studied in [TJ] we also find that 
the equations for the one-body and pair correlation functions of the system are closed. (In 
fact there are closed equations for each order of the correlation.) We define the vector p , 
Pi =< Xi >, I = 1, 2 . . . 2N and the pair correlation matrix 



u z 



at 1 - 

Z V 



(4) 



where the N x N matrices u, z and v are given by ui m = (qiq m ) , vim = (piPm) and z\ m = 
{<liPm) • It follows then from Eq. Q p and c satisfy the following equations of motion: 

dp f dp\ 

dc A „ _ r - / dc\ 

— = — ac — ca + a + — , (5) 

where the last terms in the above two equations arise from the flip dynamics and are given 
by: 

dp 
d> > 




and vd is a diagonal matrix with matrix elements [vd]u = Vu = (pf)- 

In the steady state, dp/dt = which implies p = 0. Setting dc/dt = gives the following 
set of equations for the pair correlations in the NESS: 

z T = -MzM- 1 , 

v = ^(Ml + $wM) + ^{MzTM- 1 + M- l fz T M) , 

(Mm$ - MM) + (MzTM ■ - M~ x tz T M) + 2\(Mz - z T M) = , 

($z + | T $) + (M^f?) + wfM- 1 ) + 4A(£ - = 2TT . (7) 

Using the fact that u and v are symmetric matrices we have A^ 2 + iV(iV + 1) unknown 
variables and there are that many independent equations above. 

Now consider the case of heat conduction across a harmonic chain with Hamiltonian given 
by (II]) and self consistent reservoirs attached to all sites. This is in addition to the two end 
reservoirs at fixed temperatures Tl and Tr. Each of the side reservoirs is a Langevin bath 
with a friction constant 7^ and a temperature T/, I = 1, 2, ...N, which is self-consistently fixed 
by the condition that there is no net flow of energy into the reservoir [IS] . The stochastic 
equations of motion of this system are: 

^ = -$i m qm - —Pi - ^Pi + (2 7 T i ) 1 / 2 r ?1 (t) + {2 1 ' 1 T[fl\ l {t) 
dt mi m i 

^ = -$i m q m ~ ^Pi + {2i l T' l )y%{t) I = 2, ...,N- 1 , 
dt mi 

^ = -$Nmq m ~ —Pn ~ ^Pn + (2 1 T R ) 1 / 2 r ]N (t) + {2i N T' N ) 1 l\ N {t) , (8) 
at rriN m N 

where r]i,r] N and Q, I = 1,2, iV are independent Gaussian white noise sources with unit 

variance. It is immediately established that the probability distribution P(x) in the NESS 

of this model is a Gaussian. The self consistency condition for zero current into the side 

reservoirs is given by T[ = vu = (pf) /mi. Making the identification 7^ = 2Xmi, it is seen 

that the equations for the pair-correlations in the steady state corresponding to the above 

equations are given precisely by Eq. ([7]). 

The self-consistent model was first studied by Bolsterli et al [IB] who introduced the self- 
consistent reservoirs as a simple scattering mechanism mimicking anharmonicity and which 
might ensure local equilibration and the validity of Fourier's law. The model was later 
solved exactly by Bonetto et al [18] who proved approach to local equilibrium and validity 
of Fourier's law for the ordered case, i.e where all the m^s are equal. They also obtained 



an explicit expression for the thermal conductivity of the system in all dimensions, (see Eq. 
(pb below) 



III. BOUNDS ON GREEN-KUBO CONDUCTIVITY 

Bernardin [3] considered a model of a disordered harmonic chain with a stochastic noise 
that changes the momentum of neighboring particles while keeping the sum of their kinetic 
energies constant. He obtained an exact result for the Green-Kubo conductivity of an 
ordered chain and also rigorous upper and lower bounds for the conductivity of disordered 
chains. Here we use Bernardin's approach for our model to obtain an exact expression for 
the ordered chain. We also obtain bounds for the conductivity of the disordered chain which 
are slightly improved from those of Bernardin's. 

The time evolution of the phase space density is given by Eq. ^ which we rewrite here 
in a more abstract form for convenience. 

where L = A + XS 



f-f mi oqi ^ dpi 

1=1 m=l 



n 

SP(x) = J2[P(-,-Pi,-)-P(-,Pi,-)} • (9) 
i 

The total current which is carried entirely by the Hamiltonian part can be written in the 
following form: 

k N 

J=n'Y]—(.Qi+i-Qi-i) >9o = QN ,Qn+i = 9i. (10) 
i^i mi 

The Green-Kubo expression for the thermal conductivity at temperature T is given by: 

i r°° 

k gk = lim lim — / dte~ zt (J(0)J(t)) 

1 

>6 n^o NT 2 



lim lim / dt e~ zt \ dx J e Lt ( J P eq 



o 



lim lim -^=(J, {z-L)~ x J). (11) 

z~ 5>0 N^oo iV 1 



where we have used the notation (/, g) = J dxf(x)g(x)P eq for any two functions /, g of phase 
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space variables x = (q±, . . . , qN,Pi, ■ ■ ■ ,Pn) and P eq ~ exp[—/3H] where H is given by the 
periodicized version of Eq. (Jl]) with k' set equal to 0. 

We note the following relations which are easy to prove: 

AJP eq = V ^(q l+1 - qi ^)P eq 
' mi 

i,j 

and SJP eq = -2JP eq (12) 



A. Green-Kubo conductivity for equal mass ordered case 



For the equal mass case Eq. (12) gives AJP eq = 0. This is true with or without pinning 
and corresponds to the fact that for periodic boundary conditions the current operator 
commutes with the Hamiltonian. Hence we get: 

1 f,^ 1 



z- 



Kgk = hni lim rn9 Ar / dx J — ^—^ J P eq = Jim 



^on^ooT 2 N z + 2\ y jv^oo2AT 2 A^ 



Using the form of J in Eq. (10) we then get 



kD . . 

KGK = 77T — , (13) 



where D 



4k 



2k + k + [{k ){4k + K)} 1 ^ • 

This agrees with the result of [18] [Eq. (4.18)] where the conductivity was defined as 
k = lim A r^ 00 (J iV )iV/(T L - T R ) where (J N ) = k(qipi +1 /m l+1 ) (I = 1,2,...,N - 1) is the 
average heat flux in the NESS of the system with self-consistent reservoirs specified by Q. 
The mapping between the noisy dynamics model and the self-consistent reservoirs model 
(with the transformation 7 = 2Am ) implies that for our noisy model also k = kqk in 
the ordered case. Following the methods in [18J it is easy to show that the value of k is 
independent of boundary conditions for the ordered case and while not proven we expect 
this to be true also for the disordered case for N — > 00 at fixed A > 0. In fact there is 
every reason to believe that whenever the Green-Kubo formula for kqk converges to a fi- 
nite value when iV — > 00 then it will agree with the conductivity in the NESS defined as 
k = limAr^oo lim Ti _> Tfl ^ T N(J N )/(T L - T R ). 
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B. Upper and lower bounds on the Green-Kubo conductivity 



We now consider the random case where the masses are independently chosen from some 
distribution. Bernardin's proof that the conductivity Kqk is bounded away from zero and 
infinity is based on an identity between (J7", (z — V)~ X J) and a variational expression (Eq 
(15) in 0), 



(J, {z - L)- X J) = Sup{2(u, J)-{(z- XS)u, u)-{{z- XS^Au, Au)} 



(14) 



where the supremum is carried out over the set of smooth functions u(q,p). The derivation 
of this formula is straightforward for A=0. More generally we can consider a symmetric 
L, e.g one corresponding to a stochastic dynamics satisfying detailed balance with respect 
to P eq . Then we have that u is the solution of the equation Lu = J , and both sides of 



Eq. ( 14 ) are equal to (J7", u) . For a derivation of Eq. ( 14 ) in the case L = S + A with A 
antisymmetric, see [22J. 

Lower bound: Choose a test function u = nY2iPi(<li+i ~ Qi-i)i where \i is a variational 
parameter: 



((z-XS)u,u) = (z + 2X)(u 2 ) = (z + 2X)fi 2 Tj2™i((<H+i-<li-i) 2 ) ■ 



where D is defined in Eq. (13). Denoting by [...] an average over disorder we then get: 

H 2 D[m] 



[((z-XS)u,u)\ = NT 2 (z + 2X)'- 



k 



Similarly, 



t(z - XS^Au, Au) = (z + 4A)~ 1 ((Au) 2 ) = (z + AXr'^T 2 V ( —) 

*Y \mi m l+1 J 

and averaging over disorder gives 



mim+i 





' 1 " 










m 





Thus we have: 



NT 2 



[{J,{z-L)- l J)\ > Dfi-Cfi 2 



2XD[m] 1 , r . 
where C = 1 — r{[Tn\ 

K ZA 



1 

m 



I) . 



(15) 
(16) 
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The minimum of the bound occurs at /i = D/(2C) and this gives: 



[*gk] > ^ • (17) 



Upper bound: By neglecting the last term in Eq. (14) which is clearly negative, we get 
the upper bound: 

(J,{z — L)~ l J) < ( z + 2A)- 1 ( l 7 2 ) = (^ + 2A)- 1 T^-((g m -^ 1 ) 2 ) • (18) 
Hence, 

1 kD 1 

M = M*U J >( Z - L) ~ lj) ^8\ [ m ] ' (19) 



Combining (17) and (19) gives: 

kD r . kD r 1 , / \ 

SAHCl + fc ^H-VH ) ^ ^ ^ 8A [ m ] (2 ° } 

As A — > oo, both bounds behave as 1/A while for A — > 0, the upper bound diverges while the 
lower bound goes to linearly in A. The behavior of kqk and of k in the NESS when A — > 
after A — > oo is thus not determined by these bounds and remains an open problem for 
both the pinned and unpinned random mass case. What we do know is that, if A — > with 
N finite then there is a significant difference between the pinned and unpinned cases [TJ [2] . 
As already noted in the introduction, for the pinned case all phonon modes are localized 
with a fixed localization length independent of iV and the current decays exponentially with 
system size. In the unpinned case the low frequency modes are extended and the current 
has a power law decay with an exponent that depends on the boundary conditions used pQ , 
k n ~ A" 1 / 2 for fixed BCs US] and as k n ~ N 1 / 2 for free BCs [H |2Q]. With the addition 
of the noisy dynamics which conserves energy but not momentum we expect as noted earlier 
that the conductivity k will be equal to kqk and thus strictly positive for any A > [21J. 
In the following section we evaluate (Jjv) as a function of A and N numerically and via 
computer simulations, to obtain information about its behavior when A — > 0. 



IV. RESULTS FROM NUMERICS AND SIMULATIONS 

We study the dependence of the heat current in the NESS on the system size and on 
the strengths of the disorder and noise. In all our computations we set k — 1. The masses 
{mi} are chosen from a uniform distribution between 1 — A to 1 + A. This gives [m] = 

9 



1, [1/m] = l/(2A)ln[(l + A)/(l — A)]. The average heat current from site I to I + 1 is 
given by = k(qiPi+i/mi+i) ■ In the steady state this is independent of I and we denote 
ji+1,1 — (Jn)- We note that ( Jjy) = kzij+i/mi+i and hence we can obtain accurate numerical 
values for the current in the disordered system by solving the equations for the correlation 
matrix i.e Eqs. (|7|). This involves solving large dimensional linear matrix equations and we 
have been able to do this for system sizes less than A = 512. For larger sizes we performed 
nonequilibrium simulations and obtained the steady state current by a time average. For 
small sizes we verified that both methods agreed to very high accuracy. The number for 
disorder realizations was 100 for A < 64, and varied between 2 — 16 for larger sizes. The error 
bars in our data presented below are calculated using the results from different realizations. 

A. Pinned case 

This corresponds to the case with k Q > and here we also set k' = 1. All results in 
this section were obtained by numerical solution of Eqs. ([7j. In Fig. we plot (Jn)N/ AT 
versus A for different values of the flipping rate A and with A = 0.8, k D = 4. In all cases we 
see a rapid convergence to a system-size independent value which then gives the conductivity 
k of the system. In Fig. (|2j) we plot k, obtained from the large- A data in Fig. 0, as a 
function of A. For comparison we also plot the lower and upper bounds for the Green-Kubo 
conductivity given by Eqs. 
can be thought of as a balance between the flips delocalizing the phonons and acting as 
scatterers of phonons. 

In accord with the bounds we find that at large A, k ~ 1/A while at small A, the numerical 
results suggest k ~ A. We note that for A = 0, all phonon modes are exponentially localized 
within length-scales i ~ (fc A 2 ) -1 . One can argue that for small values of A there is diffusion 
of energy between these localized states with a diffusion constant ~ £ 2 X. This leads to the 
k ~ (k 2 A A )~ l \ and we now test this numerically. In Figs, 
which suggests the scalings k ~ A~ 4 and k ~ k a ~ 2 ' 5 , which are roughly consistent with the 
expected behavior. The reason for the discrepancy could be that we are not yet in the strong 
localization regime where the prediction is expected to be most accurate. 



17fl9). It is seen that k has a maximum around A ~ 0.5. This 



3|4) we show the numerical data 
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B. Unpinned case 



As noted above for the unpinned case with A = 0, the two different boundary con- 
ditions(BC) namely fixed BCs with k' > and free BCs with k' — give respectively 
(Jn) ~ N~ 3 / 2 [3] and (Jn) ~ N~ 1 ^ 2 [4J. The difference in the asymptotic behavior of the 
current for different BCs can be understood as arising from the dependence on BCs of the 
transmission of the low frequency modes which carry the current pp. For any A > however 
we expect that the system should have a unique finite value of the conductivity, indepen- 
dent of boundary conditions, the same as for kqk- Physically we can argue as follows: The 



unpinned system without disorder has a finite positive k given by Eq. (13), which is hide 



pendent of BCs [see comments at end of sec. (Ill A)]. The low frequency modes are weakly 
affected by disorder hence we expect that as far as these modes are concerned the unpinned 
system with and without disorder will behave similarly. Since these are the modes which 
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FIG. 1: (Color online)Plot of JN/{Tl — Tr) versus N for for different values of A. The parameter 



values were set at k n 



4,k = 1,A = 0.8. All the data shown here were obtained from exact 



numerical computation. 
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FIG. 2: (Color online) Plot of K versus A obtained from the numerical data in Fig. Q. The lower 
and upper bounds for kqx given by Eqs. (17|19) are shown by the dashed lines. 
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FIG. 4: (Color Online)Plot of k versus k for different values of A and with A = 0.8 and k = 1. 
We also show a straight line with slope —2.5. 

led to the dependence on BCs for the case A = 0, we expect that for A > they will not 
have any effect. 

We now present results of our numerical and simulational studies of the unpinned chain 
with free and fixed BCs. The numerical results are obtained by solving Eqs. ([7]). The 
simulation involves evolving the system with the Hamiltonian part, the momentum flips at 
all sites and the Langevin baths at the boundary sites. For N < 512, the numerical method 
was employed to arrive at the solution for the NESS whereas for larger values of N, we 
performed simulations to obtain the data. In Fig. ^ we plot JN/ (Tl — Tr) versus iV for 
different values of A for both fixed and free BCs. The disorder strength is A = 0.8. For 
both BCs we can see flattening of the curves at large system sizes for the parameter values 
A = 0.1,0.01 implying a finite k, which is independent of BCs. For A = 0.001,0.0001 it 
appears that reaching the asymptotic limit requires larger system sizes. Using the large- N 
data in Fig. ^ we estimate the conductivity k = JN/ AT and this is plotted in Fig. (|6j). In 
Fig. ([7]) we show a typical plot of the temperature profile for the case with fixed BC (this 
was obtained using exact numerics). The profile is close to linear consistent with the fact 
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FIG. 5: (Color Online)Plot of JN/(Tl — Tr) versus N for the unpinned case with both fixed 
(dashed lines) and free BCs (solid lines) for different values of A and parameter values k = 1 
and A = 0.8. The data for N < 512 were obtained using exact numerics and in all these cases 
simulations give very good agreement with the numerics. For ./V > 512, the data were obtained 
from simulations alone. 

that the conductivity is temperature-independent. We do not see any significant boundary 
temperature jumps since the system size is sufficiently large. In Fig ^ we have the profile 
for N=128 for the case A — > which shows considerable jump in the temperatures across 
neighboring sites. It appears likely that for all A > the conductivity k is independent of 
BCs. However this is difficult to verify from simulations since one needs to study very large 
system sizes to reach the correct asymptotic limit. The reason for this can be roughly seen as 
follows. In the ordered case the conductivity k ~ 1/A and this can be understood in terms of 
an effective mean free path £ ~ 1/A for the ballistic phonons because of scattering from the 
stochastic process. Hence we can expect that, to see diffusive behavior for the low frequency 
ballistic modes, important in the disordered case, requires one to study sizes N ~ £ or 
N ~ 1/A. Finally we observe from Fig. (6) that at small A, the conductivity appears to be 
diverging as 1/A 1 / 2 . In the absence of noise the localization length £i ~ 1/w 2 , hence it is 
expected that all modes with £ L < £ or to > A 1 / 2 stay localized. The low frequency modes 
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FIG. 6: (Color Online)Plot of k versus A for the unpinned system obtained from the numerical 
data in Fig. The lower and upper bounds for kqk given by Eqs. ( 17|19 ) are shown by the 
dashed lines. Also shown is a straight line with slope —1/2. 



TABLE I: Values of Correlation Functions for n=4 (Tl = 4, Tr = 1) 



Correlation 


A = 0.1 


A = 2 


A = 10 


<p\ >-3<p\ > 2 
<p'f> 2 


0.006 


0.014 


~ io- 4 


<p*>-3<p 2 > 2 

<p\> 2 

<pj>-3<pj> 2 

<Pz> 2 
<pj>-3<pl> 2 

<p'i> 2 


0.054 


0.058 


0.049 


0.061 


0.110 


0.098 


0.038 


0.018 


0.002 


<phi>-<ii><pi> 


0.002 


0.004 


~ 10~ 5 


<q 2 Xp 2 > 


<p 2 q 2 >-<q 2 Xp 2 > 


0.011 


0.014 


0.013 


<q\><p\> 


<phl>-<il><p%> 

<q\><p\> 


0.012 


0.025 


0.023 


<p\<i\>-<ql><p\> 

<q'l><p'l> 


0.012 


0.034 


~ io- 4 
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FIG. 7: (Color Online)Plot of temperature profile (Tj = (pj/nii)) for the unpinned case (with mass 
disorder A = 0.8) with fixed BCs for N = 128, A = 5, Tl = 4, Tr = 1. The expected linear profile 
is also shown. The data was obtained from exact numerical computation. 

< to < A 1 / 2 become diffusive with mean free paths ~ 1/A thus resulting in a conductivity 
k ~ A 1 / 2 (l/A) ~ 1/A 1 / 2 , which explains the observed behavior. 

V. DISCUSSION 

We have shown here that the stationary one-body and pair correlations in the velocity flip 
model are the same (after setting j' t = 2Xmi) as in the harmonic chain with self-consistent 
reservoirs. We have also numerically investigated the dependence of the thermal conductivity 
k of the disordered harmonic chain on the velocity flip rate A. For A — > our results suggest 
k ~ A for the pinned system and k ~ A -1 / 2 for the unpinned system. Establishing these 
results conclusively requires further work. 

We note that while for the self-consistent reservoirs model the NESS is exactly Gaussian, 
this is not so for our noisy model. Instead it will be in general a superposition of Gaussians. 
Computer simulations however indicate that the single particle distributions are very 
close to a single Gaussian while the joint distribution of xi and p\ or of pi and pj , 
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FIG. 8: Same parameters as Fig. [7] except the velocity flip A — > 0. The data shown is from an 
exact numerical computation with A = 10~ 9 . We have verified that this is close to the temperature 
profile for A = 10~ 7 and expect that it is converging to the A = value. 

j 7^ /, are essentially uncorrelated. In addition, the set of equations for the four variable 
correlation functions was derived and for small values of n, they were solved numerically. 
The exact values obtained from these numerics were found to be in close match with the 
simulation results. For n=4, and Tl — 4, Tr — 1, the results from the numerical solution 
are shown in table I corresponding to 3 different values of A. It was also observed that 
these normalized forms of correlations (that go to zero when the temperatures of the 
two reservoirs are equal or when there is no velocity nipping) have a limit when either 
the difference in temperature (Tl — Tr) or strength of stochastic noise (A) goes to infin- 
ity. When we let A — > oo, we observed that the correlations involving pi and p4 went to zero. 

We are currently investigating the O(N) corrections to the pair correlations in the 
noisy NESS . These are known to behave like 1/N for certain diffusive lattice systems 
and to contribute terms of O(N) beyond those obtained from the local equilibrium to the 
variance of the particle number in the NESS. Results of this kind are also known partially 
for the continuum case with a different kind of noise, i.e instead of velocity reversals pairs 
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of nearest neighbor particles diffuse on the circle pf + pf +1 = C. 
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